Coherent Quantum Optical Control with Subwavelength Resolution 
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We suggest a new method for quantum optical control with nanoscale resolution. Our method 
allows for coherent far-field manipulation of individual quantum systems with spatial selectivity that 
is not limited by the wavelength of radiation and can, in principle, approach a few nanometers. The 
selectivity is enabled by the nonlinear atomic response, under the conditions of Electromagnetically 
Induced Transparency, to a control beam with intensity vanishing at a certain location. Practical 
performance of this technique and its potential applications to quantum information science with 
cold atoms, ions, and solid-state qubits are discussed. 
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Coherent optical fields provide a powerful tool for co- 
herent manipulation of a wide variety of quantum sys- 
tems. Examples range from optical pumping, cooling, 
and quantum control of isolated atoms [l|, [2J and ions 
Q to manipulation of individual electronic and nuclear 
spins in solid state [4, 5]. However, diffraction sets a frac- 
tion of the optical wavelength A as the fundamental limit 
to the size of the focal spot of light Q. This prohibits 
high-fidelity addressing of individual identical atoms if 
they are separated by a distance of order A or less. In 
this Letter, we propose a method for coherent optical 
far-field manipulation of quantum systems with resolu- 
tion that is not limited by the wavelength of radiation 
and can, in principle, approach a few nanometers. 

Our method for coherent sub- wavelength manipulation 
is based on the nonlinear atomic response produced by 
so-called dark resonances The main idea can be un- 
derstood using the three-state model atom shown in Fig. 
HJa). Consider two such atoms, atom 1 and atom 2, po- 
sitioned along the x-axis at x = and x = d, respec- 
tively, as shown in Fig. QJb) . Assume that they are pre- 
pared in the ground state \g) and then illuminated by the 
probe field with wavelength A and Rabi frequency f2. For 
dcA, one cannot focus the probe on atom 1 without af- 
fecting atom 2 and other neighboring atoms. Let us sup- 
pose that f2 is uniform over the distance d. In addition, 
prior to turning on the probe, we turn on a two-photon- 
resonant spatially varying control field (e.g. a standing 
wave) of wavelength A' = 2it/k' that vanishes at x — 
(i.e. has a node) and has Rabi frequency fl c (x) » flok'x 
for k'x <C 1. Destructive interference of excitation path- 
ways from \g) and \r) up to |e) ensures that the so-called 
dark state |dark(x)) = (n c (x)\g) - Vt\r)) / ^/ tt 2 c {x) + Q 2 
is decoupled from both optical fields [7(. It is the sharp 
nonlinear dependence of |dark(a;)) on fl c (x) that allows 
for sub- wavelength addressability. In particular, for atom 
1 at x = 0, dark(x)) = —\r), so that atom 1 prepared in 
state \g) responds to the probe light in the usual way. 
On the other hand, for all x such that fl c (x) 3> O, 
dark(ir)) rj \g). The space interval around x — 0, in 



which the ground state \g) is not dark, therefore, has 
width ~ f2/(f2ofc') and can thus be made arbitrarily small 
by increasing the overall intensity of the control (oc fig). 
In particular, atom 2 at x = d prepared in \g) will not 
respond to the probe provided ^> Q/(k'd). 

This selective sub-wavelength addressability can be 
used in a variety of ways. For example, one can accom- 
plish selective state manipulation of proximally spaced 
qubits via spatially selective stimulated Raman transi- 
tions. In combination with dipole-dipole interactions, 
our technique can be used, for d <C A, to generate an 
efficient two-qubit gate between pairs of atoms. One can 
implement selective fluorescence detection Q of the inter- 
nal state of an atom if \g) — \e) corresponds to a cycling 
transition (this is possible either if \r) is above |e) or 
if spontaneous emission from e) into |r) is much slower 
than into \g)). Finally, one can perform spatially selective 
optical pumping of individual atoms. Addressability with 
(!cA will be important for arrays of quantum dots [5[ or 
optically active defects [3] in solid state, where d <C A is 
often needed to achieve coupling [§]. Moreover, our tech- 
nique enables highly desirable high-fidelity micron-scale 
manipulation of atoms in optical lattices with d = A/2 1] 
and ions in linear Paul traps with d < 5 /im 

[3] (for ions, 

small d is desirable as it accompanies large vibrational 
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FIG. 1: (color online) (a) 3-level atom prepared in state \g) 
and coupled at two-photon resonance to a spatially uniform 
probe held with Rabi frequency Q and a spatially varying 
control held with Rabi frequency Q c (x). (b) Schematic of the 
setup: atom 1, at a node of the control field, responds to the 
probe, while atom 2, a distance d away, is subject to a large 
control field f2 c (d) S> and does not respond to the probe. 
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frequencies Q). Below, we analyze in detail selective 
coherent state manipulation and then estimate manipu- 
lation errors using realistic experimental parameters. 

Before proceeding, we note important prior work. Our 
approach is an extension of incoherent nonlinear tech- 
niques used in atom lithography 0] and biological imag- 
ing ■ The nonlinear saturation of EIT response that 
forms the basis of the present work has already been used 
for the realization of stationary pulses of light ll| and 
has been suggested for achieving subwavelength localiza- 
tion of an atom in a standing wave ((l^.[l3| and references 
therein). Finally, alternative approaches to solving the 
addressiblity problem exist that use Bessel probe beams 
with nodes on all but one atom 14|, place atoms into 
traps separated by more than A 11511 . and resolve closely 
spaced atoms spectroscopically [16| by applying spatially 
varying magnetic fields [17[ or light shifts [18| ■ 

As a specific example, we now analyze in detail a 
spatially selective single-qubit phase gate, |0) — » |0), 
1) — > e l< ^|l), on a qubit encoded in stable atomic states 
|0) and |1) of one atom in the presence of a proximal 
neighbor (Fig. [5]). Consider atoms 1 and 2 that have a 
tripod configuration shown in Fig.[2ja). We assume that 
the optical transitions from the metastable states |0), |1), 
and |r) up to |e) are separately addressable via polariza- 
tion or frequency selectivity. By turning on a probe field 
with Rabi frequency ~ tt, wavelength A = 2it /h, and de- 
tuning A ^ for a time t oc A/fl 2 , we would like to 
apply a 7r-phase on state 1 1) of qubit 1 via the ac Stark 
effect. To minimize errors discussed below, we turn f2 on 
and off not abruptly but adiabatically (e.g. a linear ramp 
up from zero immediately followed by a linear ramp down 
to zero). To shut off the phase shift on the nearby qubit 
2, we apply, at two-photon resonance with f2, a spatially 
varying control field with Rabi frequency £l c (x) w Qok'x 
for k'x <C 1. We assume the control is turned on before 
and turned off after the probe pulse. For the moment we 
treat atoms as point particles and return to the question 
of finite extent of the atomic wave packets below. 

The gate error on atom 1 due to spontaneous emission 
can be estimated as rjp e ~ T7(0/A) 2 ~ 7/A, where 
Pi is the population of state \i) and where we assume 
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FIG. 2: (color online) Single-qubit phase gate on atom 1. (a) 
Atom 1 (S1 C (0) = 0) or atom 2 (£l c (d) / 0). (b) Atom 2 using 
basis states {\D), \B)} in place of {|1), |r)}. (c) Schematic of 
imperfect localization of atom 1: parabolic trapping potential 
mw 2 x 2 /2 with three lowest energy levels indicated, ground 
state wavepacket of width ao, and control field fl c (x) ~ Qyk'x. 



for simplicity that all transitions are radiatively broad- 
ened and that the decay rate of |e) is 2j. To investi- 
gate the effect on atom 2, we define dark and bright 
states for this atom as |D) = (f2 c |l) — f2|r))/f2 and 
\B) = (0|1> + Q c \r))/Q, where & = ^/9? c + ft 2 and 
n c = Q c (x = d) (see Fig. HJb)). Since \D) = |1) at 
the beginning and at the end of the probe pulse (i.e. 
when fi = 0), the phase gate will be turned off if atom 
2 remains in a superposition of |0) and \D) without any 
phase accumulation on \D) or population loss into \B). 
This will be the case provided the probe field is turned on 
and off adiabatically as compared with \B) — \D) energy 
splitting, which is equal to the Stark shift A5 = fr/A of 
\B). In the limit il c ^> f2, which we will assume from now 
on, the non-adiabatic coupling between \D) and LB) has 
an effective Rabi frequency Q ma ~ H/ (Til c ) [l9( giv- 
ing population loss from the dark state into the bright 
state of order pb ~ (Qna/As) 2 ~ (0/^ c ) 6 and hence 
an error of the same order. The errors due to the Stark 
shift Vl 2 NA /As of \D) and due to spontaneous emission 
are smaller than (0/T2 c ) 6 and 7/A, respectively. 

In the simplest case, these are the dominant sources of 
error, so that the total error is 



p e ~ (7/ A) + (n/nrf 



(i) 



Plugging in fl 2 ~ A/t and minimizing with respect to A 
gives A ~ (7T 3 «6)i/4 and p e _ [ 7 /(t-02)] 3/4 , which can 
be made arbitrarily small by increasing control intensity. 

However, other sources of error exist. For d <C A, 
dipole-dipole interactions and cooperative decay effects 
may become important [20] . Cooperative decay will not 
qualitatively change the errors since the desired evolu- 
tion is close to unitary. Assuming that we have only two 
atoms and that d <C A, taking the axis of quantization to 
coincide with the a;-axis, the dipole-dipole Hamiltonian 
can be written as Hdd = (fit • M2 — 3(/ii ■ x){p2 ■ x))/d?, 
where pi is the electric dipole operator of the ith atom. 
Since most of the population will stay in |0) and |1), 
the dipole-dipole interactions involving state \r) can be 
ignored. Then, provided |0) — |e) and |1) — |e) have 
different polarizations or sufficient frequency difference, 
H dd w - 5o (|Oe)(eO| + |e0)(0e|) - 9i(|le)(el| + |el)(le|), 
where \a/3) denotes a two-atom state with atom 1 in \a) 
and atom 2 in |/3) and where go and g\ are proportional to 
g = 7/(fco?) 3 with proportionality constants that depend 
on the polarizations of the transitions. Then a pertur- 
bative calculation shows that dipole-dipole interactions 
introduce an error ~ (fig / (fl c A)) 4 [2ll ]. 

Additional errors are associated with imperfections 
in the control field node and with finite localization of 
atoms. If atom 1 was perfectly localized at a single point, 
a residual control field at the node (f2 c (0) 7^ 0) would re- 
sult in population (f2 c (0)/f2) 2 in the dark state \D) (now 
defined for atom 1). However, even if f2 c (0) = 0, atom 1 
can still interact with the control field due to finite extent 
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TABLE I: Error budget for the single-qubit phase gate. 



ao of its wave- function. Assuming S1 C (0) < flok'aa Q, 
the error due to finite atomic extent (discussed below) 
will dominate over (O c (0)/r2) 2 . 

To analyze the problem of localization for atoms in 
optical lattices and ions in linear Paul traps, we assume 
that atom 1 sits in the ground state of a harmonic os- 
cillator potential with frequency to and, therefore, has 
spread ao = \/ h/(2mu>), where m is the mass of the 
atom, as shown schematically in Fig. |2(c). We assume 
fl c (%) — Slok'x — Q, ca (a) + a), where f2 ca = f2 c (ao) and 
a is the oscillator annihilation operator. Q c {x) will then 
couple \e, n) and |r, m) only when n = mil, where \a, n) 
denotes atom 1 in internal state |a) in nth harmonic level. 
The dominant error can be estimated by keeping only 
states 1 1 , 0) , |e,0), and \r, 1). A perturbative calculation 
shows that the two limits, in which the error is small arc: 
(a) fast limit lot < 1, in which case P e ~ (f2 ca /f2) 2 ; (b) 
adiabatic limit lot 3> 1, (f2 ca /f2) 2 , in which case a small 
change in the Stark shift of |1, 0) can be compensated by 
slightly adjusting t to yield P e ~ (Q ca /Q) 2 /(tlo) 4 . 

For atom 2 centered at x = d, we have £l c (x) = £l k'd+ 
£l c k'(x — d), i.e. the desired coupling il c within each har- 
monic level is accompanied by coupling of strength ~ Q ca 
between different harmonic levels. Numerical simulations 
show that provided fl ca < 0. 1 O c (which will always hold 
below), this coupling has an insignificant effect. 

The error budget for the single-qubit phase gate is sum- 
marized in Table |TJ In general, for a given set of experi- 
mental parameters, using fl 2 ~ A/r to eliminate f2 in fa- 
vor of A, one has to write the total error as the sum of the 
errors in Table [J and minimize it with respect to flo and 
A (we assume £Iq/2tt < 1 GHz). We will illustrate this 
procedure for three systems: ions, solid-state qubits, and 
neutral atoms. Since ion and neutral atom examples will 
have d ~ A, we take Q c = Slo for them, while for solid- 
state qubits, we take f2 c = Q k'd. We take H. ca = Qok'a , 
except for neutral atoms, as discussed below. We note 
that stimulated Raman transitions [3J] , resulting in qubit 
rotations, can also be treated in exactly the same way, 
yielding similar error probabilities. Moreover, this error 
analysis is readily extendable to spatially selective qubit 
measurements and optical pumping, as well as to dipole- 
dipole two-qubit gates for qubits separated by d <C A. 

Several approaches to control field node creation exist. 



One or two standing waves can be used to generate planes 
or lines, respectively, of zero field with field amplitudes 
scaling linearly near the zeros. If one has a regular array 
of atoms (e.g. in an optical lattice), arrays of zeros can be 
chosen to have spacing incommensurate or commensurate 
with atomic spacing, allowing to address single or multi- 
ple atoms, respectively. One can also create control field 
nodes using holographic techniques [HI , which allow one 
to generate single optical vortices (such as in a Laguerre- 
Gaussian beam) or an arbitrary diffraction-limited two- 
dimensional array of them. For simplicity, we consider 
the case when atoms are sensitive only to one polariza- 
tion of the control field (e.g. if a magnetic field is applied 
to remove degeneracies). Then the quality of a standing 
wave node in this polarization component is determined 
by the interference contrast, which is limited by the mis- 
match between the amplitudes of this component in the 
two interfering waves. On the other hand, in an optical 
vortex, if the phase of the desired polarization compo- 
nent picks up a nonzero multiple of 2ir around a closed 
loop, for topological reasons this loop must enclose a line 
(in three dimensions) where the amplitude of this polar- 
ization component exactly vanishes (see e.g. 25]). Fur- 
thermore, the Rabi frequency in an optical vortex rises 
radially from the center as |f2 c (x)| ~ Clo(x/w) , where 
w > A' is the beam waist and the topological charge I 
is a positive integer. Therefore, in some cases, the use 
of vortices with I > 1 instead of standing waves or I = 1 
vortices can improve the resolution by decreasing the un- 
desired coupling of the control to atom 1. We will use an 
I = 2 vortex for the neutral-atom example, in which case 
we take £l ca — ^lo(k'ao) 2 in error #2 in Table U 

We first analyze ions in linear Paul traps. We consider 
for concreteness 40 Ca+ '26] with |0) = \4Si/ 2 , m = 1/2), 
|1) = \AS l/2 ,m = -1/2), \e) = |4P 1/2) m = 1/2), and 
\r) = |3L> 3/2 ,77i = 3/2). Note that A = 397 nm and 
A' = 866 nm are far enough apart to ignore off-resonant 
cross coupling. Then, for r = 1 /is, to /2tt = 10 MHz, 
and d = 1 — 3 fim, errors #1 and #4 from Table Q] form 
the dominant balance, so that Eq. |T]) applies and P e ~ 
[7/(t^2)] 3/4 , which is - 1(T 4 for n /2n = 1 GHz (with 
optimal A/(2tt) - 200 GHz and f2/(27r) - 200 MHz). 
This and the next two error estimates are significantly 
lower than the errors produced by naive probe focusing. 

For solid-state qubits (e.g. Nitrogen- Vacancy color cen- 
ters in diamond [2JJ), we take ao = 0.5 nm, A = A' = 700 
nm, 7/2-7T = 5 MHz, and r = 1 fis, which, for d be- 
tween 100 nm and 20 nm, makes errors #2 and #4 form 
the dominant balance, so that P e ~ (ao/d) 3 / 2 is between 
5 x 10~ 4 and 5 x 10~ 3 . For d < 10 nm, n /27r = 1 GHz is 
insufficient to suppress the dipole-dipole error (error ^5 
in Table [!}, and the gate fidelity sharply drops. 

To analyze atoms in optical lattices, we consider 87 Rb 
with |0) = \5S 1/2 ,F = 2,m F = 2), |1) = \5S 1/2 ,F = 
l,m F = l), \e) = \hP lj2 ,F = 2,m F = 2), and \r) = \AD). 
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|4D) decays with rate 2~f r — 1/(90 ns); so to reduce the 
error ~ p r j r T ~ (Q/Qc) 2 Jri~ on atom 2 (error #6 in Ta- 
ble|T|, we choose short t = 10 ns. For o;/27r = 50 kHz and 
r2 /27r = 1 GHz, errors #2 and #6 form the dominant 

1/2 

balance, so that P e ~ (f2 ca /f2 c ) (r7 r ) ' ~ 0.01. This 
error can be further reduced by tightening the traps for 
the duration of the gate either by increasing the power 
of or by decreasing the detuning of the lattice beams. 

Our selective addressability technique has several ad- 
vantages that may enable it to outperform alternative 
all-optical addressability proposals based on the gradient 
method [lij ]. First, the nonlinear response provided by 
the dark states may potentially provide our method with 
superior error scaling. Second, in the gradient method, 
the control field typically couples states that are pop- 
ulated at some point during the gate. In contrast, in 
our method, the control field is small (ideally, vanish- 
ing) on the atom that is being manipulated, while on the 
neighboring atoms the population of level \r) (coupled by 
the control to level |e)) is negligible and becomes smaller 
as the control power grows. As a result, in contrast to 
the gradient method, our method (1) avoids unwanted 
forces on atoms due to Stark shift gradients [and hence 
prevents unwanted entanglement of external and internal 
degrees of freedom] and (2) avoids excessive spontaneous 
emission, which may take place if the control field mixes 
populated stable states with short-lived excited states. 

We now outline some new avenues opened by the coher- 
ent selective addressability technique. Although we dis- 
cussed in detail only the application of this technique to 
selective phase gates (equivalently, Raman transitions), 
it has obvious generalizations to geometric gates [28[ , flu- 
orescence detection, and optical pumping/shelving, as 
well as to the generation (in combination with dipole- 
dipole interactions and assuming d <C A) of entangling 
gates between atoms. In addition to the applications to 
atoms in optical lattices, to ions in linear Paul traps, 
and to solid-state qubits, our technique may also allow 
for single-atom addressability in recently proposed sub- 
wavelength optical lattices [29j. Moreover, a combination 
of similar ideas involving dark states and the nonlinear 
atomic response can itself be used for creating deep sub- 
wavelength-separated traps and flat-bottom traps. Fi- 
nally, better optimization (e.g. using optimal control the- 
ory to shape laser pulses) can further reduce the errors. 
Therefore, we expect this technique to be of great value 
for fields ranging from quantum computation and quan- 
tum simulation to coherent control, all of which can ben- 
efit from high-fidelity addressability at d < A. 

We thank D.E. Chang, A. Peng, J. Gillen, T. Calarco, 
S. Foiling, J.E. Thomas, and M.R. Dennis for discus- 
sions. This work was supported by the NSF, Harvard- 
MIT CUA, Packard Foundation, and AFOSR MURI. 
P.Z. acknowledges support by the Austrian Science Foun- 
dation and the EU. 

Note added: after completing this work, we became 



aware of related proposals 3(| 31, 32j to use dark state 
position dependence to achieve sub-wavelength resolu- 
tion. 
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